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I. INTRODUCTION 

Decoherence is one of the main obstacles 
for building useful quantum devices. Under- 
standing the mechanisms of decoherence and 
achieving long decoherence times is crucial for 
many fields of science and applications includ- 
ing quantum computation and quantum infor- 
mation [l[, protein dynamics |2|-(4j, dynamics 
of excitons and charge separation in biological 
complexes 04^] % and the new and rapidly grow- 
ing fields of NMR and MRI with ultra small 
(microtesla) magnetic fields [Iol - il2j . In the lat- 
ter case, the Larmor frequencies of the spin pre- 
cession become relatively small (in the kHz re- 
gion) , causing the effects of 1// noise become so 
important that noise suppression must be used. 

In many situations the influence of noise can 
be modeled by an ensemble of two- level systems 
or fluctuators [lljl - iilSj . Depending on the distri- 
bution of parameters of the fluctuators, such 
as amplitudes and switching rates, and cou- 
pling constants, this model can describe both 
Gaussian and non-Gaussian effects of noise 
1 1 511 . Recent experiments with Josephson qubits 
[20|-|22(, on the quantum dynamics of excitons 
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in light-harvesting antennas in photosynthetic 
complexes demonstrated these important 
contributions of noise and thermal fluctuations 
to decoherence, relaxation processes, and quan- 
tum coherence effects. 

In this paper, we study relaxation and de- 
pha sing processes using a spin-fluctuator model 
[laUJI- In the spin-fluctuator model, fluctua- 
tions are described by a random telegraph pro- 
cess (RTP) produced by N fluctuators. Each 
fluctuator is characterized by two parameters: 
its amplitude and switching rate. Depending 
on the distribution function of fluctuators over 
amplitudes and switching rates, the RTP can 
describe noise for a broad range of frequencies 
using spectral characteristics that include both 
low- and high-frequencies noise. 

We consider the noisy environment produced 
by a large number of fluctuators, N 1. In the 
limit N — > oo, we obtain an effective random 
process (ERP) described by a continuous distri- 
bution of fluctuators. We derive a closed sys- 
tem of integro-differential equations for func- 
tions averaged over the ERP. Even though this 
system of equations is closed, it is still very 
complicated for direct analysis and even for nu- 
merical solutions. 

We study two approximations in which these 
equations are reduced to a system of differential 
equations. The first one we call the Gaussian 
approximation, because, as we demonstrate, in 
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the simplest case of a two-level system (qubit) 
under the influence of an ERP, it yields the rela- 
tion: (exp(i(p)) = exp(— (ip 2 /2)), where <p is the 
random angle of the Bloch vector. This Gaus- 
sian approximation is widely used in theoretical 
and experimental research to descre the influ- 
ence of noise on quantum systems [H, El, HH - 
[13]. In many situations, this approximation is 
very useful because (i) it captures some impor- 
tant properties of noisy dynamics and (ii) it is 
simple to apply. However, the Gaussian ap- 
proximation does not describe different "non- 
Gaussian" effects which can play a significant 
role. 

Our second approximation is based on two ef- 
fective fluctuators which include both low and 
high frequency noisy components. We show 
that this approximation goes beyond the Gaus- 
sian approach and better describes the experi- 
mental results for a superconducting flux qubit 
in a noisy environment [231 ] . 

Our main results 

• We create a new model based on an ef- 
fective random process (ERP) that in- 
cludes both slow (low-frequencies) and 
fast (high-frequencies) fluctuators. This 
model can describe the influence of noise 
on a quantum system over a wide fre- 
quency range. 

• For the functions averaged over the ERP, 
we obtain an integro-differential master 
equation which we reduce to a closed sys- 
tem of differential equations in two ap- 
proximations: (i) a Gaussian approxima- 
tion and (ii) an approximation of two- 
effective fluctuators. Both of these ap- 
proximations describe, to some extent, 
the contributions from low and high fre- 
quency noise. 

• We demonstrate that the two-effective 
fluctuator approximation accurately 
models "non-Gaussian" effects observed 
in experiments with superconducting flux 
qubits [H] . 

• We show that the two-effective- 
fluctuators model better describes the 
suppression of 1/f noise in experiments 
involving echo decay in superconducting 
flux qubits (23[. 

This paper is organized as follows. In Sec. 
II, a general model of noise based on ERP is 



introduced that describes both low-frequency 
(1//) and high-frequency noise. In Sec. Ill, we 
use the reduced density matrix approach to de- 
scribe the interaction of a quantum system with 
its environment by a master equation. For two 
cases (i) the Gaussian approximation and (ii) 
the two effective- fluctuator model, we reduce 
the system of integro-differential equations to a 
closed system of differential equations. In Sec. 
IV, the general method developed in Sees. II 
and III is applied to describe the decoherence 
of a superconducting flux qubit for free induc- 
tion decay and for echo decay experiments. In 
the same Sec. IV, we compare our theoreti- 
cal predictions with available experimental data 
and demonstrate a good agreement with experi- 
ments. We conclude in Sec. V with a discussion 
of our results. In the Appendices we present 
some technical details. 



II. DESCRIPTION OF NOISE USING A 
RANDOM TELEGRAPH PROCESS 

To describe noise we use the spin-fluctuator 
model developed in [l3|, EH- In this model, 
noise is described by a sum of N uncorrelated 

fluctuators, £n = Si=iCi(*), where Q(t) is a 
random telegraph process (RTP). The variable, 
(i(t), takes the values, — or a^. Consequently, 
Cf(*) = af = const. 




&(*)) = 0, (1) 
UtMt'))=Si j a 2 i e- 2 ^ t - t \ (2) 

The amplitude, a^, together with the switching 
rate, 7j, completely characterize the i-th fluctu- 
ator. The correlation function related to £/v(t) 
is defined as, XN(\t—t'\) = (6v(i)6v (£'))• Using 
Eqs. (P) and @, we obtain 

N 

XN{\t-t'\)=Y^a 2 l e' 2 ^ t - t \ (3) 
i=i 

Further, assuming N 3> 1, we consider contin- 
uous distributions of amplitudes and switching 
rates. The corresponding correlation function, 
x(t) = linijv^oo Xat(^), can be written as 

X(\t-t'\) = JJ dw^ye-W-'K (4) 

where, a 2 = limpj^oo Na 2 , and dw(a,j) de- 
pends on the specific distributions of am- 
plitudes and switching rates. The random 
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process described by the function, = 
lim.jv->oo ^jv(i), we call an effective random pro- 
cess (ERP). 

In order to model the characteristic be- 
havior of the spectral density of noise in 
different frequency domains, we introduce a 
family of random variables and distributions, 
{£n(X)) dw„(a, 7)}. In particular, n = 1 corre- 
sponds to low-frequency (I//) noise, and n — 
2 corresponds to the Lorentzian spectrum for 
high frequencies. (See Appendix C for details.) 

Accordingly, we introduce the ERP as £ (t) = 
J2n£ n (t)i where each £ n {t) is an independent 
source of noise. This implies (£m(*)£n(*')) = 
(to 7^ n). As shown in the Appendix A, the 
corresponding spectral density S n (uj) behaves 
as S n (uj) l/uj n in some region of frequen- 
cies. The total correlation function is a sum of 
the partial correlation functions, — t'\) = 
T, n Xn{\t - t'\), where 

Xn(\t-t'\)= [[dwnia^ye- 2 ^-^. (5) 



In this paper, we adopt the simple model in- 
troduced in [13[ for uncorrelated a and 7. We 
define the distribution function as 

dw n (a,j) = S(a- a n )V n (-y)dadj, (6) 

where a n is a typical value of the amplitude and 



Pn(7)d7 = 46(7c„ - 7)0(7 ~ 7, 



d7 



(7) 



Here, Q(x), is a step-function; and r y m , n and 
7 Cji are the lower and upper switching rates, 
respectively. The normalization constant, A n , 
is: 

I 



A n 



111(701/77^)' 
(l-7^ 1 /7?- 1 )' 



(8) 



n + 1 



In the following, we restrict ourselves to two 
important cases: n — 1 and n — 2, which are 
related to 1 // noise and to high-frequency noise 
with the corresponding spectral densities. (The 
case for arbitrary n is analyzed in Appendix 
C.) We denote j m = Tmu 7c = 7ci, and 70 = 
7 C2 . For the distribution functions, "Pi (7) and 
7^2(7), we impose conditions at the point 7 = 
7c so that 7 m2 = 7 C1 = 7 C (7™ < 7 C < 70). 

Using Eq. (|C6I) (Appendix C), we obtain, 

Xi(r) = a 1 2 A 1 (£; 1 (27 m r) - £i(2 7c t)), (9) 



X2(r) 



<J 2 A 2 



E 2 (2 1c t) £ 2 (2 7 ot) 



To 



(10) 



Computation of the spectral density, 



S n (u>) = Xn(T)COs(ujT)dT 7 (11) 





yields 



Si(u) 



5 2 (w) 



a\A\ ( ( uj \ / oj 

arctan I — arctan I 

TTLU \ V27 m / V27c 

(12) 



In 



a%An (1+uj 2 /^ 



ttuj 2 \1 + oj 2 / 



(13) 



where At = 1/ ]n(j c /~/ m ) and A 2 = 7 C /(1 
7c/7o)- 




FIG. 1: Spectral density of noise (2*y m — Is , 
27 c = l^is" 1 , 270 = 10 /xs -1 , 0-1=0-2 = 1, 
uj = 2-7t/). Red line: The total spectral density, 
S(f) = Si(f) + &(/)• Green line: Contribution of 
the slow fluctuators described by Si(f). Blue line: 
Contribution of the fast fluctuators described by 
82(f)- Blue dashed line: The Lorentzian spectrum, 
Sl(/) ~ 1/f 2 - Black line: The spectral density of 
1/f noise, S 1/f (f)=A/(2nf). 

From Eqs. (TJ2J and (JTSJ) it follows that, in 
the interval, j m < uj < j ci the spectral density 
Si(uj) describes 1/f noise. Indeed, in this inter- 
val Si(uj) « A/u>, where A = a\/(2 \n(^ c /^ m )), 
and 7 m and 7 C are related to the infrared, 
uj m = 27 m , and the ultraviolet, uj c = 2j c , fre- 
quency cutoffs, respectively (Appendix C). For 
52 (w), we obtain the following asymptotic be- 
havior: S 2 (uj) ~ 1/uj 2 (uj ^> uj c ). Thus, asymp- 
totically S 2 (uj) has a Lorentzian spectrum. 

To estimate the relative contributions of dif- 
ferent processes for low- and high-frequency 
noise, we evaluate the relation, S 2 (lu) / Si(uj) at 
frequencies uj ~ and uj as ui c . A simple com- 
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putation yields the following rough estimate: 
of 7 m ln(7 c /7 m ) 



0. 



7c 



(14) 



Taking values typical for superconducting 



qubits, 2j m 
tain 



Is and 27 c ps 1 /is , we ob- 



lO-V^/df, wwO, 
10 crj/cr 2 , w PS w c 



(15) 



Thus, for low-frequency noise the main contri- 
bution near oj — is provided by slow fluctu- 
ators (SF) with the switching rates, 7, being 
in the interval (7 m ,7 c ). However, for high fre- 



quencies, ui 



> 



the contribution of fast fluc- 



tuators (FF) with 7 > j c dominates. (See Fig. 

HI) 



III. MASTER EQUATION FOR 
AVERAGED DENSITY MATRIX 

We consider a quantum system governed 
by the Hamiltonian, H(t) (generally time- 
dependent), depending on control parameters, 
Xi ( external flux, biased current, critical cur- 
rent, etc). The noise associated with fluctua- 
tions of these parameters is described by ran- 
dom functions, 5Xi(t). For simplicity, we re- 
strict ourselves to only one fluctuating param- 
eter, SX(t), denoting it as £(i). Generalization 
for many parameters is straightforward. Ex- 
panding the Hamiltonian to first order in £(i), 
we obtain 



(16) 



To include the effects of a thermal bath, we 
use the reduced density matrix approach lead- 
ing to the master equation: 

^ = - l [H(t),p}+£p, (17) 

where the superoperator, C, describes coupling 
to the bath. 

Using (fT6|). one can recast Eq. (fT?]) as 

dp(t) 



dt 



-i[Uo(t),p(t)]+Cp(t) 



7 iK(*)V(*),p(*)]._ (18) 
For the average density matrix, this yields 

d(p(t)) 



dt 



= -i[U Q {t),(p{t))]+C(p{t)) 

-i[V(t),(X(t))], (19) 



where (X(t)) — (£,(t)p(t)), and the average ( ) 
is taken over the random process describing the 
noise. 

As before, we assume that fluctuations are 
produced by the ERP, so that £(i) = J2 n £n{t), 
and the correlation function can be written as a 
sum of the partial correlation functions, ~ 
t'\)=Y,nXn{\t-t'\). (SeeEq.EJ) 

Eq. p^|) . an integro-differential equation, is 
rather complicated. However, in two important 
cases (the Gaussian approximation and the 
approximation by effective fluctuators) we 
obtain a closed system of first order differ- 
ential equations (Appendix B). Below we 
summarize our results for £(t) = £i(t) + &(t), 
where £i(£) is related to slow fluctuators 
leading to 1/f noise, and is related to 

fast fluctuators leading to high-frequency noise. 

The Gaussian approximation. Applying the 
method described in Appendix B, we find that, 
in the Gaussian approximation, the master 
equation can be recast as follows 



d(p(t)) 
dt 



-i[H (i), </>(*)>] +£</>(*)) 

-[V(t),[K(t),{p(t))}} + 0(\\ V|| 4 ), 

(20) 



where JC{t) = J* dt' X (t - t')W (t)V(t')U(t). We 
denote V(t) = V (t)V \t)U^ (t) , and 



/ i [Uo(t')dt\ 

U(t)=T(e J ° ) 



(21) 



The approximation by two effective fluctua- 
tors. In the approximation by two effective 
fluctuators, the set of slow, £i(t), and fast, 
£2(^)5 fluctuators is approximated by two effec- 
tive fluctuators: one for SF and the other for 
FF. The total correlation function, x(K — — 
Xi(\t — t'\) + X2{\t — t'\) , is approximated as 



Xn (\t-t'\)^a* n z e 



1,2, (22) 



where a* and 7* (the effective amplitude and 
switching rate) are defined as follows: a* 2 = 
Xn(0) and 7* = -(l/2)d]n X (t)Mlt=o. (For 
details, see Appendix B.) 

Applying the method developed in Appendix 
B for an arbitrary system of stochastic first- 
order ordinary differential equations, we obtain 
from Eq. p^|) the following closed system of 
ordinary differential equations: 
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-( P (t)) = -i[H (t), (p)] + C(p) - i[V(t), (Xi(t))] - i[V(t), (X 2 (t))}, (23) 

j t (X!(t)) = -2 1 * 1 {X l (t))-i[H (t),(X 1 (t))]+C{X l (t)) 

-iaf[V(t),{p(t))}-i[V(t),{X 12 (t))}, (24) 
|<X 2 (t)) = -2 1 * 2 {X 2 {t))- l [U«{t),{X 1 {t))]+C{X 2 {t)) 

-iaf[V(t),{p(t))}-i[V(t),{X 12 (t))}, (25) 
|<*ia(f)> = -2( 7l * + ^)(X 12 (t)) i[U (t), (X 12 (t))} + C(X 12 (t)) 

-iaf[V(t), (X 2 (t)) - iaf[V(t), (X^t))], (26) 



where (X^t)) = (£1 (*)/>(*)>, (X 2 (t)) = 
(&(t)p(t)) and (X 12 (t)) = (6(*)6(*)p(t)). 



IV. NON-GAUSSIAN NOISE AND 
DECOHERENCE IN A 
SUPERCONDUCTING PHASE QUBIT 

In this section, the general method devel- 
oped in Sees. 2 and 3 is applied to describe 
relaxation effects in a superconducting qubit. 
The effective Hamiltonian for a superconduct- 
ing qubit can be written as [23[ (see also refer- 
ences therein), 

H(t) = -±fl(t) ■ v, (27) 

where er denotes the Pauli matrices. We as- 
sume that H(t) depends on the control param- 
eters, \i, of the system, including external flux, 
biased current, critical current, etc. Limiting 
ourselves to a single fluctuating parameter, A, 
and expanding the Hamiltonian in Eq. (|27[) to 
first order in the fluctuations, S\(t), we obtain, 

U{t) = -\n.a-\sX{t)^--a, (28) 

where, for simplicity, we assume that fi does 
not depend on t. In the eigenbasis of the unper- 
turbed Hamiltonian, Eq. (29) takes the form 

H(t) = - ^na z - ^D x ,z5\(t)a z 

- lD x ,xSX(t)a ± , (29) 

where D\. z = dfl/dX and a± denotes the trans- 
verse spin components, either a x or a y . (We 
adopt the notation of Ref. [H|.) 

Below, in the framework of the ERP model, 
we obtain the relaxation rates, and compare our 
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results with the results which follow from the 
well-known Bloch-Redfield (BR) theory [HHil 
applied to the external noise [18|]. Before pro- 
ceeding, we present here some important results 
of the BR approach. 

In BR theory, the dynamics of a two-level sys- 
tem is described by two rates: the longitudinal 
relaxation rate, Ti = T± , and the transverse 
relaxation rate, T 2 — TZ . BR theory is valid 
if Ti,T 2 ^> t c , where r c is the fluctuation corre- 
lation time. The transverse relaxation rate, T 2 , 
is a combination of Ti and the so-called "pure 
dephasing" rate, T v , 

T 2 = ^T 1 +T V . (30) 

In terms of the spectral density of noise, 5*a(w), 
these rates are defined as follows [l8j : 

^ =nDl ± S x (n) 7 (31) 
T v =nDl z S x (0). (32) 

In our approach, fluctuations of the parame- 
ter, A, are described by an ERP. Thus, 5X(t) = 
£n(t). Further, we restrict ourselves to 
consideration only the case, n = 1,2. Then, 
SX(t) = £i(i) + £ 2 {t), where £i(i) describes the 
contribution to the ERP of SF, and £ 2 (£) de- 
scribes the contribution of FF. The spectral 
density of noise can be written as S\ (lu) = 
Si(uj) + S 2 (uj). 

Since only FFs have small correlation times 
and satisfy the conditions of applicability of BR 
theory, we use the spectral density of FFs given 
by Eq. (Ti~3")) to calculate the relaxation and de- 
phasing rates provided by the BR theory. We 
obtain 
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T 1 =KDl ± S 2 {n)=Dl ± 



of 7c 



In 



f> 2 (l- 7 c/7o) Vl + ^747o 



1 + r» 2 /4 7c s 



r v> =wDi z S 2 (0) = Di z -± + ^ 



47c 



7o 



Note, that the validity of the BR theory is 
restricted by the condition: riT2,r 2 T2 <C 1, 
where r 2 = (l/27 c )(l +7 c /7o) is the effective 
correlation time of the FF. 

The above effective rates can also be obtained 
directly from the averaged expressions for the 
partial rates, 

J 2 

?i = D l±J ^T^dMl), (35) 

7c 

}° 2 

T V =D\ Z J ^-dw 2 { 1 ). (36) 

7c 

A. Pure decoherence 

Let us consider the Hamiltonian ([25)) for a 
pure decoherence case. Then D\,± = 0, and 
the Hamiltonian rl(t) takes the form 

U(t) = ~Sl(t)<r a - ^D x , z 6X(t)a z , (37) 

where D Xz = dSl/dX. 

The equation of motion for the density ma- 
trix, ip — \H(t), p], reduces to only one compo- 
nent: 

^■poi = iQ(t)poi + iD x , z 6X(t)p 01 . (38) 

The matrix elements, poo an d pu, are constant. 

Setting poi(t) = p i(0e wo(t) ; where ip (t) = 
J Q fldt is a regular phase, we find that poi(t) 
satisfies the following differential equation: 

-p m (t)=iD x , z 5\(t)p 01 (t). (39) 
at 

Its solution can be written as 

Poi(t) = e^*V(0), (40) 
t 

where ip(t) = J D x ^ z 5X(t)dt is a random phase, 
o 

After averaging over the random process, we 
obtain (p 01 (t)) = (e tv< -^)p Q1 {0). This yields 



(33) 
(34) 

I 

(poi(i)) = e^°( t )(e^(*)) / 9 i(0). Thus, the prob- 
lem of obtaining exact solution of Eq. (f38|) 
reduces to the computation of the generating 
functional, (e lip( ^). 

Returning to Eq. (1551 . one can see that, for 
averaged components of the density matrix, it 
takes the form 

j t iPoi(t)) = iQ(t){p 01 {t))+iD XiZ {SX(t)p 01 {t)), 

(41) 

In what follows, we obtain solutions of Eq. 
(|4I[) in the Gaussian approximation and in the 
two-effective-fluctuator approximation. We ap- 
ply these solutions to describe two widely used 
experimental protocols: (i) free induction decay 
and (ii) echo decay. We compare our theoreti- 
cal predictions with the experimental data (23j , 
and demonstrate that the experimental results 
(i) are described by the Gaussian approxima- 
tion and (ii) that the details of the dynamics 
of the signal decay can be understood by us- 
ing slow and fast effective fluctuators. We also 
demonstrate that the approach based on two 
effective fluctuators allows one to fit the exper- 
imental data better. 



1. Gaussian approximation 

In the Gaussian approximation Eq. (|4"Tj) can 
be presented in the form 

j t {poi(t)) =iQ(t)(p 01 (t)) 
t 

-Dl^J X {t-t')dt'){ Pol {t)), (42) 
o 

where x(t — t') = (5X(t)X(t')). Its solution can 
be written as 

(p Q1 (t)) = e'^Wrt^X^Wlmo)), (43) 

where (po(t) = f Q(r)dT is a regular phase, 
t 

ip(t) = D Xz J SX(t')dt' is the random phase ac- 
o 
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cumulated during the time, t, and 



evolutions, so that I13J, 111 



& 2 (t))=Dl z J J X (\t'-t"\)dt>dt", (44) 
o o 



, t/2 t 

tp(t) = D x , z ( J 5X(t')dt' - J 5\(t')dt' 

^ t/2 



is the variance of tp(t). 

Thus, in the Gaussian approximation, the 
random phase of the free-induction decay is In the Gaussian approximation, we obtain 
Gaussian distributed, and we obtain the well- 
known result for the generating functional, 



(e^«) = e 



-(i/2)<v 2 (t)> 



(47) 



(48) 



( e Mt)\ = e -(i/2)<v"(t)) < 



(45) 



where 



Using the spectral function of noise, S\(ui), one 
can rewrite (g5| as [H, QJ] 



(c^W) = exp ( - | do^Msinc 2 ^ 

— OO 

(46) 

where sine a; = sinx/x. 

In the echo experiments, the total phase, 



(V 2 (*)) = £>L / dt'dt" x (\t'-t"\) 




t/2 f 



-4D^y df y d*" X (|f -*"!)■ (49) 

t/2 

In terms of the spectral density, the echo de- 



ij)(t), is defined as difference between two free cay can be written as [l3l [l8j 
I 



OO 

\ ( t 2 r^2 f i n l \ ■ 2 wi • 2 w * 

) = exp - — Z> Ai<s / do^Msm -jsmc — 



(50) 



In Appendix C, we obtain explicit expressions 
for {ip 2 {t)) and {^ 2 {t)). 

Using the asymptotic formulae for the expo- 
nential integrals, E n (z), [32J, we find that, for 
7 TO t <C 1 (jet < 1), the free-induction decay 
produced by SF is given by 

<e*W> = exp ( - t 2 Dl z A(ln -L + 0{l)) \ , 

(51) 

where A = a 2 i /(2\n(j c /^f m )). Substituting 
uj m — 2^ m , we find that (fBTj) is exactly the 
same expression that is used in the literature 
for estimating the quasistatic contribution of 
1/f noise to the free- induction decay [18|. In 
the same limit, for the echo decay we obtain 

(e #(4) ) =exp(-t 2 Dl z A\n2), (52) 

which coincides with the corresponding formula 
obtained from Eq. (|5U| for uj m t -c 1 [35] . 



The contribution of low frequencies, wi< 1, 
in (|45[) obtained in the limit 7 m ^,7c^ <C 1, is 

(e^^exp^-^x). (53) 

This coincides with the corresponding expres- 
sion widely used in the literature [lH • 

2. Two-effective fluctuator model 

In this Section, the SF and FF introduced 
above are approximately described by two ef- 
fective fluctuators with the following correla- 
tion functions (see Appendix B), 

Xndt-ifD^a^e- 2 ^-^, (n=l,2). (54) 
where 
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the effective amplitude and switching rate being For the averaged functions: (poi(t)}, 
cr„ and 7„. Computation yields (Xi{t)) = (£i(*)/?0i(*))) (^(t)) = (6(*)Poi(*)) 

and (Xi 2 (t)) = (Cl Poi(*)), we obtain the 



7i 



(56) following closed system of first-order differential 



ln(7c/7m)' equations: 



7c ln(7o/7c ,, 7 , 
72 = -: ; • 57) 

1 - 7c/7o 



-(poi (t)> = in(t)<poi(t)> + ^, 2 ((IiW) + (X a (t))), (58) 

^<X x (t)) = -2 7 i<X x (t)) + i^t)^*)) + i£» A , 2 «X 12 (t)) + a 2 <p i(i)», (59) 
at 

^-(X 2 (t)) = -2 72 (X 1 (t)) + iO(t)(X 2 (t)) + iD A , 2 ((X 12 (i)> + a 2 2 (p 01 (t))), (60) 

^<X 12 (i)) = -2( ll + l2 )(X 12 (t))+iQ{t)(X 12 (t))+zD XtZ (a 2 2 (X 1 (t))+a 2 1 (X 2 (t))). (61) 
ore 



The solution of Eqs. 
ten as, 



II) can be writ- 



e^WSi^SaWpoiCO), (62) 
e^WiiW^aWpoxCO), (63) 
iD x>z (X 2 (t)) = e^W^^iaWpoiW, (64) 



(A)i(*)) 
iD XtZ {X t {t)) 



Dl z (X 12 (t)) 



*i(t)*2(*)P0i(0),(65) 



where ^o(*) = Jo^{t')dt'. We denote by 

(z = 1, 2) the generating functional of the RTP 



$i(t) = (expjz / dr^CT)^}), (66) 



where v 2 = D? z a 2 . The generating functional 



satisfies the second order differential equation 

M3, 



d 2 ^ 
dt 2 



' dt 



= o, 



(67) 



with the initial conditions being $j(0) = 1 and 
d$((0)/dt = 0. 

Free induction and echo decay solutions for 
a single fluctuator. In the following, we con- 
sider solutions of Eq. (|6"T)) corresponding to free 
induction signal and echo signal experiments. 
Previously Eq. ([67| was studied in [3, US |3tJ ■ 

• The solution corresponding to the de cay 
of the free induction signal is given by fla . 



Hit) 



-at 



sinh(7i//j*) 



+e 7lt cosh(7i/Zit), 



where yn = yjl- Vi 2 /~f 2 . 

• In the echo experiments, the 7r-pulse with 
duration, t\, is applied at time, r = t/2, 
to switch the two states of qubit. It is 
assumed that t\ <C t. The correspond- 
ing solution for the functional (t), with 
the initial conditions $^(0) = 1 and 
d<f> e t (0)/dt = 0, is written as [H 



*?(*) 



-at 



■ (/ii sinh(7i/Xit) 



+ cosh( 7 i/Xjt) + /i 2 - l). 



(69) 



B. Comparison with experiment 

In this section, we compare our theoreti- 
cal predictions with the experimental data ob- 
tained in 1231) and the theoretical results of the 
model [33| . The measurement of the decoher- 
ence due to 1// noise was done for the flux 
qubit described by the effective Hamiltonian 

m 



Ho 



e A 



(70) 



with the energy difference between two eigen- 
states E 01 = + A 2 . 
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The diagonalized Hamiltonian, with the fluc- 
tuations only of Eqi, can be written as 

a 

where D\ a , z = dEoi/dX a , and the term, 5X a (t), 
describes the fluctuations of A re in the Hamil- 
tonian. In the experiments j23| . the authors 
studied the decoherence due to fluctuations of 
(i) the normalized external flux, — $ cx /$o, 
where $o is the flux quantum, and (ii) the 
SQUID bias current, The contributions 

from different decoherence sources were sepa- 
rated, so that the fluctuations of and lb were 
observed independently. 

We consider two approximations: (i) the two- 
effective fluctuator approximation and (ii) the 
Gaussian approximation. 

1. Two- effective- fluctuator approximation 



O 0.5 




12 3 4 

Time(us) 



FIG. 2: Sample A from [H. Echo decay, $ e (t) = 
(e**^). The blue solid line is fit by two-fluctuator 
solution, (e"" (t) ) = $f(t)$!(t). The green solid line 
is the theoretical predictions of Ref . [33] . Red solid 
line corresponds to the Gaussian decay, (e 1 ^^) = 
e v e > _ The experimental data (not shown) are 
obtained for decoherence at value An^ = 0.0009 
(Fig. 4a, Ref. JH). 



In the approximation of two effective fluctu- 
ators, 8X(t) = £i(t) +£,2(t), and 



&(*)&(*0> = <W 



i = l,2. (72) 



Our task is to determine the fitting parame- 
ters: (vi, 71,^2,72), where vf = D\ z af, and 



7i 



7c(l -Im/lc) 

ln(7c/7m) 



72 



7 c ln(7 /7 e ) 
f - 7c/7o 



(73) 



The switching rates, j m and 7 C , are chosen 
according to the available experimental data 
for the spectral behavior of 1/ f noise, namely, 
7 m ~ I s _1 and j c ~ I /is -1 . Then, the only 
two free fitting parameters are 70 and vi . Their 
values are chosen from the best fit of theoretical 
results to the experimental data. 

To fix the value of vi, we use the experi- 
mental data for echo decay fitted to the Gaus- 
sian decay, exp(— (r^ s (A)i) 2 ) , using the rela- 
tion from [H, HI 



r£ B (A) = 



dE Q1 



OX 



(74) 



The constant, A \, is determined from the ex- 
perimental data describing the behavior of the 
spectral density of I// noise, S\(uj) — A\/u), 
at the frequency, / = 1Hz [23j |. 

Inserting A x = o-f/(2\n(j c /j m )) into Eq. 
(Tflll. we obtain 



21n(7 c /7 m ) 



In 2 



(75) 



In Fig. P2J we compare our theoretical predic- 
tions with the experimental data obtained for 
decoherence of a flux qubit with fluctuations 
of the external normalized flux, (sample A 
from 23]). To fit our theoretical results to the 
experimental curves, we set two cutoffs for I// 
noise as: j c = 0.5 /is -1 and -f m — 0.5 s _1 . 
Then, calculating the switching rate, 71 , we 
obtain 71 = 0.04 /is -1 . From Fig. 4c (in |23l|h 
describing echo dephasing rate T 9 E vs An^, we 
extract T 9 ipE — 0.8 /is -1 , and, then, using (f75|) . 

we obtain Vi = ^Jv^ = 4.92 /is -1 . The param- 
eters, V2 and 70, are chosen by best fitting our 
curve to the experimental data. For the high- 
frequency noise we obtain the upper cutoff as 

70 = 4.25 /is" 1 . For the switching rate 72 this 
yields: 72 = 1.2 /ts _1 . The amplitude v 2 we 
choose as v 2 — 2.72 /is -1 . In Fig. P51 we com- 
pare our theoretical predictions with the exper- 
imental data obtained for decoherence in a flux 
qubit with fluctuations of the external flux, 
(sample B from (23|). The fitting parameters 
obtained in the same way as for sample A are: 

71 = 0.04 /is -1 , V! = 21 /is -1 , 72 = 5.75 /is -1 , 
v 2 = 12.45 /is -1 and T 9 ipE = 3.75 /is -1 . In Figs. 
PJ] and [3j the echo decay of the two fluctuator 
model (blue curves) resulted from both low and 
high frequency fluctuators. 

In order to determine the contribution from 
only 1/f noise to this echo decay, we present 
in Fig. |3]the decay (for samples A and B) pro- 
vided by only a slow effective fluctuator with 
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u.o u.o i 02 04 0( 

Time{tis) r/me(^s) 



FIG. 3: Sample B from [H]. Echo decay, $ e (t) = 
(e^CO), The blue solid line is fit by the two fiuctu- 
ator solution. The red solid line corresponds to 

the Gaussian decay e f> E . The experimen- 
tal data are obtained for decoherence for value 
An = -0.0008 (Fig. 4d, Ref. 




1 2 3 4 

Time(^is) 



FIG. 4: Suppression of 1// noise in echo-decay ex- 
periment in a flux qubit at fluctuations of the exter- 
nal flux. The blue solid line corresponds to sample 
A, and the red solid line corresponds to sample B 
from Ref. [11 . 



the same parameters, 71 and v±, as those indi- 
cated in Figs. [2] and [3J One can see from Fig. 
2] that in both cases, a suppression of 1// noise 
is up to 95% in the time-interval, (0 — 1/is). 

In Fig. [5J we compare our theoretical results 
for echo decay with the experimental data ob- 
tained for decoherence in a flux qubit for fluc- 
tuations of SQUID bias currents If, (sample A 
from [Hj]). In all considered cases, we find that 
our solutions based on two effective fluctuators 
better fit the experimental data than the theo- 
retical description of the Gaussian approxima- 
tion used in |23j . 

We consider also free induction decay, and 
compare the obtained effective decoherence 
rate, T^ Fl with the experimental data and the- 
oretical results of the Gaussian model [H, [23[ . 
For free induction decay, the two-fluctuator so- 



FIG. 5: Echo decay, $ e (t) = (e^ (t) ). The 
blue solid line is fit by the solution of the two- 
effective-fiuctuators model, with the choice of 71 = 
0.04MHz, vi = 10.5MHz, 72 = 2MHz and v 2 = 
50MHz. The red solid line corresponds to the ex- 
ponential decay, e~ r * Et with T^, E = 1.7MHz. The 
experimental data are obtained for decoherence at 
value of SQUID bias current It = — 0.7/iA (sample 
A, Fig. 3c, Ref. [H). 
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FIG. 6: Sample A from [23l ]. Free induction signal 
decay, & (t) = (e l,p{t) ). Blue solid line fits to the 
two fluctuators solution, $ / (t) = &{(t)&l(t). Red 
solid line corresponds to the Gaussian decay, $^ = 
e -( r ^F*) 2 w ith r» F = 3.97 ^is" 1 . Green dashed 
line presents the Gaussian approximation for free 
decay, {e iv{t) ) = e - {1/2)(ip2{ - t)) . (Data are taken 
from Ref. [2^] for decoherence at various flux biases 
n^, sample A.) 



lution is 

(e *(t)) = ${(t)$/(t), (76) 

where $f (t) (i = 1,2) are given by Eq. (I55|) . 

Expanding (1761) into a Taylor series, we ob- 
tain 

{e i ^)=l-\{vl + vi)t 2 + .... (77) 
Then, comparing with the Gaussian decay, 
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-( r ^F*) we obtain 



1 ipF 



(78) 



Substituting vi = 4.92MHz and v 2 = 2.72MHz 
(sample A), we find T 9 vF = 3.97MHz. 



Computation for the sample A of the the ra- 
tio, r^/r^, yields T^/T^ a 4.96. This is 
in a good agreement with the theoretical pre- 
diction, T 9 vF /T 9 vE < 5 [H, and with the es- 
timate from the experimental data yielding a 
ratio between 4.5 and 7.5 12311. 



2. Gaussian approximation 

In the Gaussian approximation, free induc- 
tion signal decay is described by 



/ e Mt)\ = e - 



(l/2)< V 2 (t)) 



(79) 



where & 2 (t)) = (p?(t)) + (^(t)), and (<^(i)) 
(n = 1, 2) is given by Eq. (jCfS)) . 

In Fig. [BJ we compare the Gaussian approx- 
imation (green dashed line), the two-fluctuator 
solution (blue solid line) and the Gaussian de- 
cay (red solid line). As one can see, the Gaus- 
sian approximation and the Gaussian decay 
yield practically the same results. However, the 
two-fluctuator solution shows non-Gaussian os- 
cillatory behavior. 

We also considered the echo decay signal for 
the Sample A. Our theoretical results for echo 
decay follow from (gSJ and (|39]l 




FIG. 7: Sample A [23|]. Echo signal decay, 
$ e (i) = (e^C). The blue solid line is fit by the 
Gaussian approximation, (e^ (t) > = e -(i/2)W 2 (t)> . 
Green solid line: two fluctuator solution, (e 1 ^^) = 
<1>5 (t)<t>| (*) ■ The red solid line corresponds to Gaus- 
sian decay, (e^ (t) > = e~ (r ^') 2 . 




Time(ns) 



FIG. 8: Suppression of 1// noise in echo-decay ex- 
periment. The blue solid line corresponds to the 
two fluctuator solution with only slow fluctuators. 
The red solid line corresponds to the Gaussian ap- 
proximation with slow fluctuators. 



( e ^W) =e -(i/2)<v 2 (*)> ; (80 ) 

where (^(t)) = + (^)), and (^(t)) 

(n=l,2) is given by Eq. (fC33]l 

In Fig. [3 we present theoretical results (sam- 
ple A) for echo decay for: the two-effective- 
fiuctuator model (green curve); the Gaussian 
approximation (blue curve), and Gaussian de- 
cay used in [19] (red curve). 

In Fig. [H we compare the theoretical results 
for suppression of 1/f noise by slow fluctua- 
tors in the two-effective-fluctuator model and 
the Gaussian approximation. One can see that, 
up to 1 /is, both descriptions give similar re- 
sults. For times larger the 1 microsecond, the 
Gaussian approximation does not describe the 
suppression of 1// noise, since fluctuators with 
7 » 7 C dominate. 



V. CONCLUSIONS 

The approach based on modeling of noisy 
environment by an ensemble of two-level sys- 
tems (fluctuators) is widel y u sed for quantum 
solid-state systems [rj - fl^ Il8l [l9j. Recent ex- 
periments with Josephson phase qubits [20I - I22I ] 
demonstrated the importance of noise of all fre- 
quencies in decoherence processes and stimu- 
lated theoretical discussions on the contribu- 
tions of low- and high-frequencies fluctuators 




In this paper, we have discussed the SF 
model for continuous distribution of fluctua- 
tors to describe both, low- and high-frequency 
noise. We considered two approximations of 
our model: the Gaussian approximation and 
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two-fluctuator approximation, and compared 
our theoretical predictions with the experimen- 
tal results for decoherence of a superconducting 
flux qubit [23|. We showed a good agreement 
between our theoretical model and experimen- 
tal results. 

We should emphasize that the two-fluctuator 
approximation leads to the non-Gaussian be- 
havior in the signal decay. The non- Gaussian 
effects, yielding contribution to a particular be- 
havior of the tail in the spin echo signal, are 
very strong for free induction signal decay. The 
main problem is that available experimental 
results on superconducting qubits (including 
those reported in Ref. |23J) may not have a 
good enough precision to distinguish Gaussian 
and non-Gaussian behavior. However, it is no 
doubt that the non-Gaussian behavior is rele- 
vant to many situations and can help to under- 
stand better the nature of noise and its action 
on the system under consideration. 
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Appendix A: Some properties of random processes 



1. Random telegraph process 



In this section, we derive some useful formulae for the random telegraph process (RTP) defined 
by £jv(i) = J2i=i Ci(*) w ith the correlation function, XN(\t — t'\) = (£jv(i)6v(t'))> gi yen by 

N 

Xiv(|«-t / |) = E a " e " 27<1 ^' 1 - ( A1 ) 

i=l 

We assume that the RTP is described by N uncorrelated fluctuators, Q(t). Each fluctuator switches 
randomly between the values —1 and 1 with the probability 1/2, so that (f(t) = a 2 — const, and 
after averaging over the initial states of each fluctuator, the following correlation relations hold 

[Mill 

m)) = 0, (A2) 
(C i (t)Q(t'))=S ij a 2 i e- 2 ^- t '\t>t', (A3) 
Af^ii.ia, . . .,t n ) = a 2 e- 2 ^-^M^ 2 (t 3 , . . . ,t„), (A4) 

where 

Aij(*i,t2,...,t„) = <Ci(*i).-.Ci(*n))» *i >*2 > ••• >*„. (A5) 
From Eqs. (IA2|) - (IA4|) . a recursive formula follows 

M^(t u t 2 ,...,t n )=XN(ti-t 2 )M^_ 2 {t 3 ,...,t n ), (A6) 

where 

M£(t 1 ,t 2 ,...,t n ) = (Z N (t 1 )...Z N (t n )) ) h>t 2 >--->t n . (A7) 
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The RTP is conveniently described by the generating functional (3lj . 

t 

$iv[t;w(r)] = ^cxpji I dT£ N (T)v(T)}y (A8) 



Applying Eq. (|A6I) and using the Taylor expansion of Eq. (|A8[) . we obtain an exact integral 
equation for the generating functional $Ar[i;w]: 

t ti 

$ N [t;v(T)] = l- f dti J dt 2 XN{ti-t 2 )v{ti)v(t 2 )$ N [h;v(T)}. (A9) 



One can transform this integral equation into the integro-differential equation, 

t 

j t ^ N [t-v{r)] =-v(t) J dtiXJv(*-*i)«(*i)^[*i;«(r)] (A10) 
o 

Let R[t; £n(t)] be an arbitrary functional. Then, using Eq. (|A6I) and a Taylor expansion in 
£,n{t), one can show that the following correlation splitting formula holds: 

(6v(M6v(*2)i?[i;6v(r)]) =Xff(h-t 2 ){R[t;^ N (r)]), h > t 2 > t. (All) 

To calculate the correlator (^(t)R[t; £zv(t)]) for t < t we use the following relations [3l|: 

t 

(^(t)ie[t;^(r)+»?(T)]> = (6v(i)ex P { J dr£ N (t) -^r}) R[t; t,(t)], (A12) 

o 

where ^(r) is a deterministic function. With the help of Eq. (|A10[) . we obtain 

t h 
(^N(t)R[t'Mr)+r}{r)]) = J dt lXN (t - *i)(^y ex P { ^ dr^(r^r ^(r)] 

o 

dt lXN (t - i^^-^Rfa^r) + ^(r)0(ii - r)]). (A13) 
o 

Taking the limit 77 — > 0, we find 

t 

((N(t)R[t;( N (T)]) = J ^ixjv(t-ti)(^r%ti;^(r)]), (A14) 


where 





t 



R[t,t 1 ;£ N (T)}=R[t;£ N (T)Q(t 1 -T + 0)}. (A15) 
By differentiating (|A15[) with respect to time t, we obtain 



(Ai6) 



This formula generalizes the differential formula [29-31 

d 



(_ + 2 7 )<C(W;C(t)]> =(C(*)|i2[*5C(r)]), (A17) 
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taking place for the RTP described by £(t) with switching rate, 7. 

Theorem 1. For the random telegraph process, £/v(i), the following relation holds: 

(6y(*W;6rfr)]) = ^^T^ famfrM-r)]), t'>t, (A18) 
where R[t; £n(t)] is an arbitrary functional. 
Proof. Writing £^ (f ) as 

iV iV 

JV 

we can employ the fact that C, 2 (i) = const [23,130]. Next, using the relation Xiv(0) = X) £ 4 2 Wj we 

»=i 

obtain 



&(*) 



1 W 

7^E^(*)0(*) = 1. (A20) 



WiOi xa(0) .. 
Inserting (|A20I) into the l.h.s. of Eq. (|A18[) . we find 

W)#[*;6v(r)]) = (teV)^UrmtN(T)]) ^E^(*)0(*)eiv(t')^;^(T)]). (A21) 

Then, applying (|A6I) . we obtain 

(ZN(t')R[t;Z N (T)]) = XN{t '~ t) (tN(t)R[tMT)]) - ^—J2mKj(t))(^(t')R[t^ N (r)}), H > t. 
Xiv(0) xn(0) j^. 



(A22) 



Since for, i ^ j, (Ci(*)Cj(*)) = 0, this yields 



(MOiZfofrrfr)]) = ^ ( J ] (^N(t)R\t^N(r)l t>>t. (A23) 



Corollary. In the limit A — > 00, one has 



(t( t ')R[t;ttT)}) = ^^mmaT)}, f>t. (A24) 



where x(*' - i) = lim./v_ >00 xjv(*' - *)• 



2. Effective Random Process 

We define the effective random telegraph process (ERP) for N >• 1, as £(t) = limjv^oo ^jv(t), 
considering the continuous distribution of amplitudes and switching rates. The correlation function, 
x(t) = limjv^oo Xiv(i), can be written as 

X(lt-t'l) = lim X;^ e "* Y<|t " t ' 1 = // Mmye-W-'l, (A25 ) 

where, cr 2 = limAr^oo JVo 2 , and dw(a, 7), depends on the specific distribution functions of fluctua- 
tors on the amplitudes and switching rates. The main relations for the ERP can be obtained from 
the previous section by taking the limit N 00. Below we present the most important formulae. 
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The generating functional for the ERP being defined as 

t 

$[t;u(r)] = ^expji J dT^(r)u(r)}^ (A26) 
o 

satisfies the following integral equation: 

t ti 

*[*; «(t)] =1- J dh J dt 2X (h ~ t 2 )v{t 1 )v{t 2 )^[t 2 ; w(r)]. (A27) 
o o 

One can transform Eq. (|A27[) into the integro-differential equation, 

t 

— $[t; w(t)] = -v(t) I dt lX (t ~ ti)»(ti)$[t i; «(T)] (A28) 
o 

For an arbitrary functional £(r)] the following correlation splitting formula holds: 

(((ti)((t2)R[t^(T)])= X (ti-t 2 )(R[t;((r)}}, h>t 2 >T. (A29) 
Finally, the differentiation formula (|A16I) takes the form 

t 

!<£(W;£(r)]} -y^ 1 | x (t-t 1 )^y(%ti^(r)]) = (£(i)|i?[i;£(r)]). (A30) 
o 

Relation to the Gaussian random process. We would like to mention here an important consequence 
of the central limit theorem concerning a relation between ERP and the Gaussian random process. 
Assume that for individual fluctuators the correlation relations are given by 

(&(*)}= 0, (A31) 

(Ci(t)C j (t')) = ^S ij e- 2 ^- t '\. (A32) 

Then, for N —> oo, the ERP, defined by £jv(t), becomes a Gaussian Markovian process with an 
exponential correlation function [2^, [3l[ 

mW))=^e-^ t - t '\ (A33) 

where £(f) = linijv->oo Cjv(i)- Thus, the iV-fluctuator RTP, with the same switching rates, 7, 
and the amplitudes, ct 2 /N , for a finite number, N, is an approximation of a Gaussian Markovian 
process. 



Appendix B: Stochastic differential equations 

We consider a system of first-order stochastic differential equations 

^x(t) = i(i)x(t) + £(t)B(f)x(t), x(0) = x , (Bl) 

where describes ERP, so that (£(*)£(*')) = - i') (* > t') and 

X(|t-t'|) = | dw(a, 7 )a 2 e- 2 ^-^. (B2) 

In what follows we study two approximations leading to a closed system of differential equations for 
averaged variables: (i) The effective fluctuator approximation and (ii) the Gaussian approximation. 
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1. Gaussian approximation 

In the interaction picture, we introduce the new variable x(t) = C/ _1 (i)x(i), where 

t 

/ / A(t')dt\ 

U{t)=T[el )x(t), (B3) 

with a T-ordered exponential on the r.h.s. For x(i), Eq. (|B1I) takes the form 

j±(t) = i£(t)B(t)±(t), x(0) = xo, (B4) 

where we set iB(t) = C/ _1 (<)B(t)C/(i). Eq. (|B4) can be recast as 

t 

|x(t) = J(i)£(i)x(0) - B(t) | atW)B(t% *(t')dt'. (B5) 

o 

After averaging over ERP, we obtain the following integro-differential equation 

t 

|<x(i)) = -B(t) | *(t - t')B(t')(i(t'))dt'. (B6) 
o 

For practical purposes, Eq. (IB6|) is not very useful. However for some reasonable assumptions, 
it can be simplified. First, employing (|B6[) one can write 

, t" 

(x(f )) = (*(*)) - | 4 dt'B(i') | X (f - t")B(t")(Sc(t'))dt". (B7) 
* o 

Then, inserting (|B7I) into Eq. (|B6[) we obtain 

t 

|(x(i)) = - J(t) | X (t - t')B{t')dt'{m + 0(\\ B || 4 ). (B8) 
o 

Considering the last term at the r.h.s. of Eq. (IB1|) as a perturbation, one can approximate Eq. 
(|B6[) as follows: 

t 

|(x(t)) = -B(t) | x (t - t')B(t')dt'(x(t))- (B9) 
o 

Its formal solution can be written as 

t t' 

(*(*)) = T{ exp ( - / Vjdf | x (f - t")B{t")dt") } (x(0)> 


t i 

^{ exp ( - 1 1 1 B(t')x(t' - t")B(t")dt'dt") } (x(0)) . (BIO) 



o o 



As can be seen, it has the form of solution for the Gaussian random process [29M31l | . 
Returning to (x(i)), we obtain the first-order differential equation 

±{x(t)) = A(t)(x(t)) + B(t)V(t)( X (t)), (Bll) 
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where 

V(t)= [ dt'xit-t^U^Bit^U-^t). (B12) 
Jo 

As an illustrative example, let us consider the following stochastic differential equation: 

^_x(t) = iAx(t) + iv£(t)x(t), (B13) 

with A and v being const. Its solution can be written as follows: 

(x(t)) =e ivo ^(e i ^)(x(0)), (B14) 
t 

where ifo(t) — At is the regular part, and tp(t) — v J x(t — t')dt' is the stochastic phase accumulated 

o 

at time t. 

In the Gaussian approximation, we find that the average (x(t)) satisfies the differential equation 

t 

j t (x(t)) = iA(x{t)) - « 2 ( J X (t- t')dt') (x(t)). (B15) 

o 

The solution of Eq. (|B15|) is given by 

(x(t)) = (e Uf *V>)e-"V>(e ivl % (B16) 

where 

t t' t t 

«(*) = v 2 f dt' [ X (t - t')dt" = \v 2 ( ( X (\t - t'\)dt'dt" = \{ V 2 (*)). (B17) 



2 7 7 2 

bo oo 

From here and (|B14j) . it follows that the decay law for (e lv ^) is the Gaussian, 

(e^W) =e -<v 2 ( t ))/ 2 . (B18) 
This agrees with the general conclusions made in this section. 

2. Two-effective-fluctuator approximation 



Averaging Eq. (|B1|) over the ERP, we obtain 

d 



dt x x(t)) = A(t)<x(t)) + B(t)<Xc(t)), (B19) 
where (Xf(i)) = (£(t)x(t)). Using (jBll) . and taking into account that (Xf(0)) = 0, we obtain 



t t 
(Xf(t)>= /^^i(0(X ? (i')>^'+ fxV-l')Bil'\ xt/'w//'. (B2DI 



Taking the derivative on both sides of Eq. (|B20[) , we obtain 
|(X ? (t)) =i(tXXs(t)) + x (0)B(t)<x(t)) 



1 / 9 x (t - t') 



x(o) 7 at 



(i(t')(X c (t'))' + x(0)S(t')(x(t')>)*'- (B21) 



o 
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Finally, we obtain the following closed system of integro-differential equations: 

!<x(i)) =i(i)(x(t)) + B(t)(Xs(t)), (B22) 

|<X ? (i)) =i(t)(X 5 (t)> + X (0)B(t)(x(t)) 
t 

+ ^o) / dXi m ° (Mt')(xdt')) + x(o)%'K*(0))^ (B23) 



In this section, we consider the system of Eqs. (|B22[) . (|B23I) in the approximation that the ERP 
can be approximated by a random telegraph process with the correlation function, x*(l^ ~~ 

X {\t-t'\)= [ dw(a n )a 2 e- 2 ^ t - t '^x*{\t-t'\) = a* 2 e- 2 ^ t - t '\ (B24) 



where the time-independent parameters, a* and 7* (the effective amplitude and the switching rate) 
are defined by the following expressions: a* 2 = x(0) and 7* = — (l/2)<91nx(i)/<9i| t= o. 
To proceed, consider Eq. (|B23[) rewritten as 

|<X s (t)) =A(t)(X i (t))+ x (0)B(t)( x (t)) 
t 

+ J dlnxtt - t') x{t _ t>) j-_]_ i(f0(X£( ^ )) + B( t ')(x(t>)))dt>. (B25) 


Usually lnx(i) is a slowly-changing function. Then, replacing d In 7 (t — t')/dt by its value at time, 
t = t' , one can approximate the integral on the right side of Eq. (|B25|) as follows: 

t 

^^^ Xit - t>)(-^A(t')(X ( (t>))> + B(t')(x(t')))dt> 


t 

* 91n ^~ f/) L=t Jxft-O^CO^^+^OW))^. (B26) 


Inserting (|B26|) into (|B25I) . and employing (IB20[) we obtain 

!<X s (t)) = i(t)<X ? (0) - 2 7 *<X € (t)) + a* 2 B(t)(x(t)), (B27) 

where 7* = -(l/2)<91nx(*)/<9*|t=o and a* 2 = x(0). Next, combining (jBT9| and (jB27| . instead 
of a system of integro-differential equations, we obtain a closed system of first-order differential 
equations 

|(x(i)) = i(i)(x(i)) + B(t)(X4(i)), (x(0)> = x(0), (B28) 

|<X £ (i)) + 2 7 *<X S (*)> = Mt)<X*®) + a* 2 B(t)(x(t)), (X C (0)) = 0. (B29) 

This system of differential equations describes RTP with the amplitude a* , switching rate 7* and 
the correlation function [3l[ 

X*(|*-*'|) = a*V 2 ^ (B30) 
In Fig. IH1 we compare the exact correlation functions, 

Xi (t) = a 2 A 1 (^ 1 (27„ i r) - £i(2 7c r)), (B31) 

^.W*^-*^). (B32) 
V 7c 7o / 
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with their approximated expressions, Xn ~ °n ex P(~ ^Int) given by ([51)1 . The choice of parameters, 
7 m and 7c , was motivated by the range of frequencies for 1// noise. The parameter, 70, was chosen 
to better fit both exact and approximate correlation functions. Note, that the correlation function 
in (|54p which describes the low frequency noise, Xii is n °t very sensitive to variations of the 
parameter, 7 m . Further, when fitting the experimental data, the parameters, j m and j c , were 
essentially the same as in Fig. [SJ As can be seen, the approximation (IB24|) describes the behavior 
of the exact correlation functions reasonably well for the region of parameters which we use. 




0123401234 
Time(^s) Time(ns) 



FIG. 9: Correlation functions, Xnif)-, (blue line) and exponential correlation functions, Xn(t) = exp(— 27„i), 
(red line). Upper panel: Low-frequency noise defined by Xi(*) (7m — 0-5 s _1 , j c = 0.5 /is" 1 ). There is 
good agreement up to ~ 2/is. Bottom panel: High- frequency noise defined by x^[t) (7c = 0.5 ^is" 1 , 
70 = 2 /is -1 ). In all cases a\ = 1. 

The system of Eqs. (|B28[) . (|B29|) approximately describes an ERP by RTP defined by a single 
fluctuator. Below, we will describe a model with two effective (low- and high-frequency) fluctuators. 
The advantage of this approach is that we calculate in a straightforward way the coefficients a* 
and 7*. 

Two- effective- fluctuators model 

Let us consider the same system of first-order stochastic differential equations as above 

^x(t) = i(t)x(t) + £(t)B(t)x(t), x(0) = x , (B33) 
at 

with the RTP described by two uncorrelated fluctuators, (i(t) and 2£ 2 (i), so that £(t) = C2(t)+C2(t)> 
and 

(Ci(*))=0, (B34) 
(Ci(*)Cj(0> = Sija^e- 2 ^-^, i = l,2. (B35) 

We set (Xi(t)) = (Ci(*)x(t)), (X 2 (t)) = (^f*W*)> and (X 12 (i)) = (Ci(t)Ca(*)x(t)>. Applying 
the formulae of differentiation for an RTP (29Hl3l| , we obtain the following system of differential 
equations for averaged variables: 

|<x(t)) = i(t)< X (t)) + swax!^)) + (x 2 (t))), 

|<Xx(t)) = -2 7 i<X 1 (t)> + i(*)<X 1 (*)> +B(t)«X 12 (t)> + af<x(f)», 
|<X 2 (t)> = -2 72 (X 1 (i)) + i(*)<X 2 (i)) + J3(i)«X 12 (i)) +a|<x(i))), 
^<X 12 (i)) = -2( 7l + 72 )<X 12 (t)) + A(t)(X 12 (t)) + B(t)(a 2 2 (Xi(t)) + a?<X 2 (t)». (B36) 
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Appendix C: Properties of the correlation functions 

We consider a family of random variables and distributions, {£ n (t), dw n (a, 7)}, in which each 
£ n (i) describes an independent ERP: (£ m (t)£n(t')) = (m ^ n). Then, the total correlation 
function is a sum of the partial correlation functions and x(|t — t'|) = XmXnd* — *'l)> can ^ e 
written as 

(l*-*'l)=E [ [ dw^a^e-^r (CI) 



We define the distribution function, dw n (a, 7), as 

dw„(cr, 7) = <5(cr - a n )V„(^)dad^, (C2) 
where, er n , is a some typical value of the amplitude, and 

V n (l)dj = A n Q( lcn - 7)6(7 -Tmjjp n=l,2,..., (C3) 

here, 6(a;), denotes the step-function, 7 mn and 7 Crl are the lower and upper switching rates, 
respectively. The normalization constant given by, 



1 

71=1 



ln(7c 1 /7m 1 ) 



A, = < 



TT~t-K ' 



(C4) 



I (1-7V/75T 1 ) 

is obtained from the normalization condition, J dw n (a, 7) = 1. 
Inserting (|C3|) into (jClj) . we obtain, 



Xn{\t-t'\)=alJJV n { 1 )d 1 e- 2 ^- t '\. (C5) 

From (|C5j) it follows that a 2 — x(0), and straightforward computation yields, 

( , 2 A ( E n {2 lmn T) E tl {2 1c ,t) \ 
Xn{r) = <J n A n — — x — , (C6) 

\ 7™„ 7c„ / 

where E n (z) denotes the Exponential integral [Hj]. 

It is convenient to describe each noise source by its spectral density, 

00 

S n (u) = ~ J Xn(r) cos(wr)dr, (C7) 


00 

and, as it can be easily seen, a 2 = 2 J S n (uj)du}. Employing Eqs. (|C1|) and lC7|) . one can obtain 



the following integral representation for S n (ui): 

S n {u) = - [ 2 ^ dw n (a,j), (C8) 

7T J ^1 + U 



where 



^(°) = - 71T 2 ( C9 ) 

7T 47^ + U! 2 



is the Lorentzian spectral density of the fluctuator with the amplitude a and switching rate 7 [13| . 
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Performing the integration in Eq. (|C7|) . we obtain for n > 2 

_!)fc+i ( | 



1 I(n+l)/2] 

S-„(w) = -alA n 1 n - x V 

7T ■ ■ 



fe=l 



(n - 2fc)w 2 



2fe \ in-1k 



n-2k 



1 / l + (o;/6») 2 

2 n U + (WCn) 2 



ra = 2p, 



arctan 



(CIO) 



(Gil) 



arctan I — ) , n — 2p + 1, 



where b n = 2 r y mn and c„ = 27 Cn . For n = 1, 2, the computation yields 



Si(u>) = 



U! \ ( U) 

arctan | — J — arctan — 



7YUJ 

a\A 2 f l + (u;/b 2 ) 2 

S 2 w = 2- In 

7rw^ \1 + (cj/c 2 r 



(C12) 
(C13) 



We impose on the distribution functions "Pi (7) and ^2(7) boundary conditions at the point 
7 = 7c so that 7 m2 = j Cl . Further, we denote j m = j mi , 7 C = j Cl , and 70 = 7 C2 (7m < 7c < To)- 
Using these notations, we obtain 



S 2 {u) 



arctan 



2j„ 



7YUJ 

_ ojM / l + c 2 /4 7c 2 \ 
ttw 2 ^l+ W 2/47o 2 /' 



( w 
— arctan 

V2 7c 



This yields the following asymptotic behavior of S\(lS) and S 2 (uj): 



27T7 m ln(7 c /7„ 
2w ln(7 c /7 m ) ' 

2cr 2 7 e (l - 7m/7c) 

7rw 2 ln(7 c /7 m ) 



1 _ 2hl 

7c 



u < 27 m , 

27 m < u) < 2 7c 

w > 2 7c , 



and 



47T7 C 



To 



2cr 2 7c 



tt(i - IcIloY 



u < 2 7c < 2 7o , 
In I ), 2 7c <o;< 27 , 



a; 



27c 



2 ct|7c ,/7o, ., 
— In — , u > 2 7o . 

I tt(1 - 7c/7o \7c 



(C14) 
(C15) 



(C16) 



(C17) 



5sH 



o-f 7m ln(7c/7 mj / 1 7c 
of ln(7c/7m) 



To 



1+cj 2 /4 7 2 

I a 2 2(l- 7c / 7o )- Vl+^ 2 /4 7o 2 



In 



u « 0, 



w £ 2 7c . 



(C18) 



From Eqs. (|T2"j) and (|T3|). it follows that in the interval, 7 m < w < 7 C , the spectral density S\{u)) 
describes 1// noise. Indeed, in this interval S\{w) ~ A/oj, where A — a 2 /(2ln('y c /j m )). For S 2 (uj) 
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we obtain the following asymptotic behavior: S2 (w) ~ l/u> 2 (10 ^> io c ). Thus, asymptotically S^w) 
yields the Lorentzian spectrum. 

Writing the spectral density for 1// noise as Si/f(ui) = AQ(u> c — u;)0(cj — oj m )/u>, where uj c and 
uj m are ultraviolet and infrared cutoff, respectively, we obtain 

S 1 (Lu)diu^2 f S 1/f (u)du = 2A\aiuJu m ) = g 2 ln (^M") . ( C19 ) 
Jo ln(7c/7m) 

From here it follows: "f c /jm ~ w c /uj m . Thus, 7 m and j c are related to the infrared and ultraviolet 
frequency cutoff, respectively. Further we assume u c = 2"f c and uj m — 2j m . 

As can be seen from Eq. (jClip , our model covers various asymptotic aspects of the spectral den- 
sity, S(u) = S n (u>), including 1// noise and the Lorentzian spectrum as some particular cases. 
This allows us to include into consideration the more complicated behaviors of the spectral density. 

Estimates of correlation times for superconducting qubits. Following (29l l3lj . we define the 
correlation time related to Xn{ T ) as 



1 



" Xn(0) 

From here, employing Eq. (IC6[) . we obtain 

1 - h/d 



Xn(r)dr. (C20) 



For b n <C c n , this yields 



61 ln(ci/6i) ' 

(n-l)(l-(b n /c n ) n ) 
nb n {l-{b n /c n ) n - 1 ) ' 

1 



n = 1, 

(C21) 



bi ln(ci/&i) ' 

(C22) 

n - 1 

nf 1. 



nb n 

Using Eq. (|C20|) . we calculate the correlation time of 1/f noise to be 



27™ ln(7 c /7 m ) 

For 7 m <C 7c, this yields 



n = 1 ] m/7c . (C23) 



27mln(7 c /7 m ) 
Computation of the correlation time r 2 yields 



T i « o — rr / \ • (C24) 



r 2 = -Ml + ^). (C25) 



47c V 7o 
From Eqs. (|C23|) and (|C25|) we obtain 

-<^ln(7c/7™)- (C26) 
n 27 c 

For superconducting qubits various experiments demonstrate that the frequency interval of 1// 
noise is / ~ (1Hz - 1MHz) [II]. Substituting 27™ = Is -1 and 27 c = l^s -1 into (IC24I) . we 
obtain an estimate of the effective correlation times as t\ ~ 0.01s. The experimental data on the 
ultraviolet cutoff of the spectral density are unknown, so 70 is unknown parameter. Supposing 
7o 3> 7c, one can estimate the effective correlation time as t 2 ~ l/(47 c ). Once again, assuming 
that 7c ~ 0.5/Lts -1 , we obtain r 2 ~ 0.5^s. So, the fluctuations due to &(t) have a shorter correlation 
times than fluctuations related to 1// noise, r 2 <C t%. Thus, indeed, the SF produce mainly noise 
with the spectrum ~ 1/w, and the FF lead to the spectrum ~ 1/lu 2 . 
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1. Free induction signal decay 

For a superconducting qubit in the Gaussian approximation, free induction signal decay is defined 

by (e l¥>(i) ) = e - (1/2) ^ 2(t) \ where tp(t) = D x<z J 8\{t')dt' is the random phase accumulated at time 

o 

t, and 

t t 

{^{t))=Dl z J j Xx {\t'-t"\)dt'dt". (C27) 
o o 

The correlation function, xx{ T ): 01 the ERP defined as SX(t) — En can be written as the 

sum of the partial correlation functions, X\( T ) — En XnW, and the overall accumulated random 

i 

phase, ip(t), is given by <p(t) = E„ -Da.z I €n(t')dt'. From this we obtain (ip 2 {t)) = E n (#U*))i 

o 

where 

t t 

<*£(*)} = Dx,z J J Xn(\t' - t"\)dt'dt". (C28) 



Computation of yields 

' E n+2 {b n t) E n+2 {c n t) 1/1 1 \ t ( 1 I 



(cp 2 n (t))=2 n Diy n A n \ 



z/H-1 n+1 



iCr 1 b^+O + nG" c«)) 



(C29) 



2. Echo decay 



In echo experiments, the total phase, ip(t), is defined as the difference between two free evolutions 

t/2 t 

r/,(t) = D x , z J SX(t')dt' - D x . z J 5\(t')dt'. (C30) 

t/2 

In the Gaussian approximation, one obtains (e^W) = e^ 1 / 2 )^ 2 ^, where 

, t t t/2 t v 

<V> 2 W) = £>L f 1 1 A'df'xAd*' - - 4 J dt ' J dt "^(\t' - f'l)) ■ (C3i) 

^0 t/2 ' 

Inserting XA (|f - = E„ Xn(\f - *"|) into Eq. g^TJ, we obtain (?/> 2 (i)) = E(V4W), where 

71 

(t t t/2 t v 

1 1 dfdt" Xn (|f " ~ A J dt ' J dt "Xn(\t' - t"\)l (C32) 
t/2 ' 

Computation yields 

,2 2, I .E n+2 {b n t/2) E n+2 (c n t/2) E n+2 (c n t) 



\ Cn Cn 

E n+2 (b n t) , 3^1 1 ^ , t ^ 1 1 ^ j( , T;! 
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